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Abstract 

We apply the Sudakov technique to description of the multiple 
production of lepton pairs in peripheral collisions of ultrarelativistic 
heavy ions. For heavy ions with Zia ~ <C 1 one needs a care- 
ful treatment of the multiple Coulomb exchange between colliding ions 
and screening effects, whereas interaction of real or virtual lepton pairs 
with colliding ions can be neglected. We demonstrate that while the 
inclusive spectra are modified by multiple Coulomb exchange between 
the colliding ions, the Coulomb corrections to the momentum inte- 
grated multiplicity distributions do vanish. After transformation to 
the impact parameter representation the probability of n lepton pair 
production is shown to obey the Poisson distribution. The relevant 
cross section is obtained. 



1 Introduction 

The multiplicity and the distribution of lepton pairs produced in the Coulomb 
fields of two colliding relativistic heavy ions are closely connected to the 
problem of unitarity. When heavy ions collide at relativistic velocities their 
Lorentz contracted electromagnetic fields are sufficiently intense to produce 
a large numbers of such pairs. Usually the process of lepton pairs production 
is considered as pair creation in the classic Coulomb potential of a charge 
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moving along a straight line. Such an approach allows one [|I| to investigate 
the impact parameter dependent total probability of the pair creation P{b), 
which by definition is connected with the total cross section (J = J P{h)(Ph. 
As was noticed in 0] the probability of single pair production calculated to 
lowest order in fine structure constant ol = ^ = (throughout we use 
units for which c = h = 1) aX small impact parameters exceeds one, thus 
violating unitarity. This excess begins at impact parameters smaller than 
the Compton wavelength of the electron Ac = — = 386 fm and at energies of 
practical interest (RHIC & LHC). 

Allowance for the finite size of colliding nuclei doesn't remedy the situ- 
ation, because that would affect only the impact parameters comparable to 
tOhe nuclei radii, which are much smaller than the Compton wavelength of 
the electron R <^ Ac- 

As was shown in this problem can be solved by taking into account 
the possibility of multiple pair production, whose relative contribution grows 
with energy and dominates at small impact parameters. Since this early 
publication there has been much work has in this direction (see e.g. and 
references therein) with the common for all statement: the probability to 
produce n lepton pairs in the Coulomb field of heavy ions colliding at fixed 
impact parameter h can be approximately represented as a Poisson distribu- 
tion i.e. P(6,n) = ^^(,-w{h)^ ^Yieie W(b) is the average multiplicity of 
pairs at a fixed impact parameter. 

Because of the somewhat controversial situation in the subject and its im- 
portance for the operation of relativistic heavy ion colliders (RHIC & LHC), 
in the present communication we revisit the multiple pair production based 
on the powerful Sudakov technique of calculation of high energy Feynman 
diagrams. Recently it has been applied p to the calculation of Coulomb cor- 
rections to single lepton pair production in relativistic heavy ion collisions. 
Here we extend the Sudakov approach to get the probability of multiple pair 
(n > 2) production in relativistic heavy ions collisions. For heavy ions with 
charge the numbers satisfy the conditions Zia ~ Z20. <^ 1, Z1Z2Q; ^ 1, 
one needs full allowance of multiple Coulomb interaction of colliding nuclei, 
whereas the secondary interaction of produced pairs (real or virtual) with 
the Coulomb fields of colliding ions can be neglected. 

The paper is organized as follows. In section ^ we introduce the nec- 
essary definitions and discussed briefiy the Sudakov'derivation of the Born 
amplitude for single pair creation. Then We show how one can sum the per- 
turbative series to get the amplitude for the process of n pairs production in 
a compact analytical form. 

In Section ^ we consider the Coulomb interaction between the colliding 
ions and show how it affects the amplitude and cross section of the process 



2 



Pi 



p'l 



A, 













rig S5 




) Tie 






J ^ 


















; ^ e+ \ 


; ^ e+ 













p-i 




2+ 





P2 



A 



Fig. 1: The peripheral process of creation of n lepton pairs with photon 
exchanges between nuclei Ai, A2 and light-by-light scattering blocks. 




under consideration. Apart from the conventional multiple photon exchange 
between ions, this interaction include also the so called "screening correc- 
tions" which are the result of the ions interaction through the virtual lepton 
pairs and are responsible for the probability of vacuum-to-vacuum transi- 
tion. In the language of the Abramovski-Gribov-Kancheli unitarity rules 
0, the diagrams with lepton pair loops can be associated with the QED 
model for the pomeron, and the n-pair production amplitudes can be associ- 
ated with the unitarity cuts through n exchanged pomerons. In section | we 
report the multiple production amplitudes in the impact representation in 
which the summation of the perturbation series for the n-pair production can 
be done in a compact closed form. Finally, using this amplitude we obtain 
the probability of n pair production which can be cast in the Poisson form, 
and the total cross section of all pair production process. 



2 The amplitude of the n pairs production 
process. 

The typical Feynman diagram (FD) describing the n lepton pairs production 
in the collision of relativistic nuclei with atomic numbers Ai,A2, with He 
exchanged photons between colliding nuclei as well as screening effects e.g. 
the insertions of light by light (LBL) scattering blocks is drawn in Fig. |l|. 
Upper and lower blocks in Fig. ^ describes many virtual photons interac- 
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tion amplitudes with nuclei. They contain the complete set of {ne + n + 2ns)\ 
Feynman diagrams. To avoid the multiple counting in what follows, we will 
multiply the relevant amplitude by the factor l/{ne\nj{2\)'^'^"). 

As was mentioned above we restrict ourselves to ions with charge numbers 
such that 

Zi^2 > 1, Zia ~ <t: 1, ^1^20; > 1, (1) 

which permits us to omit the multiphoton exchanges between the produced 
pairs and colliding ions. The dominant mechanism is a production of a single 
lepton pair per collision of equivalent photons, i.e., one lepton loop per two- 
photon ladder. The alternative mechanism of multiple pair production per 
collision of equivalent photons, i.e., multiple lepton loops per two-photon 
ladder, is suppressed by inverse powers of Z1Z2. 

For the description of a peripheral process of n lepton pairs creation i.e. 
the process 

Ai{Zi,pi)+A2{Z2,p2) Ai{Zi,p[)+A2{Z2,P2)+e+e4n) + - ■ ■+e+e_(r„), 

n = qX + qi (2) 

It is convenient to use the Sudakov parameterization for the four-momentum 
of all exchanged photons (for details see P) 

ki = aiP2 + PiPi + ki_i, d'^ki = -daid(3i(fki±, (3) 

S={Pl+P2), S Pi = M- :$> m , pi=pi-p2 — , P2=P2-Pl — , 

s s 



Pi=P2 = [ ^ ] , piki± = p2ki± = 0, s = 2pip2 = 2pip2. 



S2 



Here pi are light-like four- vectors build frompj, Mj are the masses of colliding 
nuclei, m and s are the electron mass and the total center mass energy. 

The denominators of intermediate states of the nucleon Green functions 
for upper and lower blocks are the same and have the following form 

($^ki)' + «0, sJ^A- (5Zki)' + ^0. (4) 

i i i i 

Peripheral process is characterized by small values of longitudinal Su- 
dakov parameters a^, (3i and the transverse momenta of the order of electron 
mass 

\ai\ ~ lAI < 1, -k1^ ~ (5) 
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Further simplification follows from the form of the nominators of the ex- 
changed photons Green function (we work in the FejTiman gauge). Using 
the Gribov's representation for the metric tensors 



9,^ = + - ( Pl'P2'' + Pl>~2^ ) . (6) 



it is easy to show that for the typical conversion of the nuclei currents 
Jfi{pi)J'^{p2) only one term, which contains the scalar products of a nucleus 
current with a four-momentum of another nucleus, becomes relevant (with 
power accuracy) 



(7) 



It can be seen that the quantity J^{pi)J^{p2)/s remains finite with the large 
values of s. This fact provides great simplification of the spinor structure of 
the amplitude 

u{v'l)P2{Vl + X1 + Mi)p2 ...{P2 + XN + Mi)p2u{pi) ^ S^+^Ni, (8) 
U{P2)P1{P2 + + M2)pi ...{P2+VN + M2)piU{p2) ^ S^+^A^2, 

Ni = ^U{p[)p2u{pi), N2 = ^U{P2)P1U{P2). 

Besides we have ^lA?"!!^ = ^|iV2p = 2 for the nuclei with the spin 1/2 and 
{Nil"^ = |A'"2p = 1 for the scalar one. Using the identity 

N 



j:^^ ^=n- (9) 

perm ^ ^ \ ^ i=l 



one can be convinced that the amplitude describing the upper and lower 
blocks in Fig. |I| can be put in the form 



h = N,\]( + I, = N2f]( + ^^1 

1=1 ^ ' 1=1 ^ ' 

(10) 

with = np + 2n, + n — 1. 
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Ml^> = is{%7^afZiZ2NiN2 Z.l ' ^ B^p = v{q+)0^pu{q_), (12) 



This expressions contain all dependence on Sudakov parameters a^, 
(the 4-momenta of exchanged photons in the peripheral kinematics in denom- 
inators of their Green functions can be considered as Euclidean two-vectors 
kf = sa,[3i + kj^ ^ k^^ = -ki^). 

At this stage the integration over Sudakov parameters can be done, be- 
cause the dependence of the amplitude on ctj, /3j provides the convergence of 
the relevant integrals 

„ Af r ^ 

hH dai = {2mfNi, hJJ d[3, = {2mfN2. (11) 

i=i i=i 

Let us now consider the single pair production. The amplitude of the 
process (1) in its lowest order (Born approximation) reads |^ 

Bai3 is the Compton tensor |^ for pair creation by two virtual photons with 
polarization vectors ei(gi), e2{q2)- Q'i(2) are the 4-momenta of exchange pho- 
tons and r = g+ -(- g_. Strictly speaking the squares of these 4-vectors qf 
do not vanish in the limit qi^ 0. This fact becomes essential when one 
calculates the total cross section of a single pair production process. For the 
case of two or more pairs production (which is our case) the replacement 
qf = — qi^ can safely be done. 
Using the gauge invariance 

q'^B^f, = q^B^/, = (13) 

one can perform the replacement 

= ^ qi q2 , = r, Si = Sa2/5l, (14) 

s Si |qi| 

The quantity Si is related to the square of the invariant mass of a pair 

si = {qi + q2? = rsi- (qi + qa)' = (g+ + q-f. (15) 

Two-dimensional vectors Cj can be interpreted as a polarization vectors of 
exchanged virtual photons. 

Using (|T^|T^) one can rewrite the Born amplitude (|T^) in a form 

M^] =isNiN2B{qi,q2), (16a) 
Biqi, q2) = {^7iafZiZ2 ^f''}f\ (16b) 

■si|qi||q2| 
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Now we are able to construct the amplitude for the process of n pairs pro- 
duction. Bearing in mind the expressions ([TI|) the matrix element of two 
pairs production can be represented as the convolution of two Born terms 
from expression ( p,6b| ) 

Ml^]=t^sN,N2 J Bik,k-n)Biq-k,q-r2-k)^. (17) 

A straightforward generalization to the matrix element in the case of n pairs 
production reads 

/n—l / J2u \ 

l[(^Bih,h-n)^jB{h,h-rn), 

n-1 

h = q-J2k,. (18) 
1=1 

Thus one can see that the amplitude for multiple pair production is solely 
determined by the convolution of the amplitudes corresponding to single pair 
production. 

In the language of the AGK unitarity rules, this result can be interpreted 
as unitarity cut through all exchanged pomerons. 



3 The Coulomb exchanges between ions 

Let us now consider the effect of m photon exchanges between nuclei Ai, 
A2. The arguments given above leads to the following matrix element for the 
process of n pairs production with m photons exchanges among the colliding 
ions 

/n— 1 / (i^k- \ n—l m 

JJ {B{ki,ki-ri)^-^ \ B{kn,kn-rn), k^ = q -^ki - ^Xi- 
i=l ^ ^ i=l i=l 

Another effect which we take into account is the possibility of the ion- 
ion interaction through the LBL blocks (screening effect in Fig. 1), which 
in the AGK language H] is equivalent to the exchange by additional un-cut 
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Fig. 2: Typical kernel of the ion-ion interaction through the LBL blocks. 



pomerons. It is associated with the iteration of a typical kernel 

L [ Y{l)d% (20a) 



iZTl'^ J |li||l2||ll — l||l2 + 1| Si 

p = n"^^'5e°(/i)e^(/i - l)ej{l2)ei{k + I), (20c) 

where we rearrange the "extra" phase volume of longitudinal Sudakov pa- 
rameters a2, Pi in terms of invariant mass square of LBL block and extract 
explicitly the boost degree of freedom of LBL block 

da2dPi f d(3i f dsi f dsi 



s{a2piy J Pi J sf J s{ 
L = In7i72, si = (/i + > 4ml 

The structure (20) entered the matrix element (19) in the compact form 

-n(^>^(ii)A) (22) 

1=1 

Real part of ( pOc[ ) (equal to one half of its s-channel discontinuity) is 
related with (|16b| ) 

g?n"^^^e?(/i)ef (/i - I)el{l2)el{l2 + ^) = \j Bapih. r - /i)e°(/i)ef (r - h) 

■ B^sih -l,r + l- h)el{li - l)el{r + I - h)d^r, (23) 
where d^r is the phase volume of the intermediate pair 

(27r)22e+2e_ ^ ^ 

We will show later that the imaginary part of P is irrelevant either for the 
total cross section or for the probability of n pairs production distribution. 
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4 The impact parameter representation for 
the ampHtude 

The last step in building the matrix element for the process of n real pairs 
creation consist in transformation of the obtained above expressions in the 
impact-parameter representation and summation over all eikonal photons 
and LBL blocks. For this we introduce the identity 

n-l 



52(^k„ - q + + ^Xi + XI lm)c?'k, 

i=l i=l 1=1 

n—1 



e-'^i^exp 



(k„ + X ki + X Xi + 5^ Im) 

i=l 1=1 i=l 



TT TT 



(25) 



Using this expression the summation in fig and Ug can be easily done with 
the result 



with 
A{P) 



'-NiN2 I e-^'»V*(''^)e-^[^^'')/2+^^(^)]]^S(p,ri; 



i=l 



TT 



, n>2, 



(26) 



.kpU^^ (27) 



+ z^{p) = J F(l)e^'^^, B{p, rO = J B{k, k-r,,. 

*(pA) = -aZ,Z2 / , , = -2aZ,Z2Ko{Xp) (28) 

J X + ^ TT 



where Kq{\p) is modified Bessel function (Mac-Donald function). 

The phase volume of final state which consists from the scattered nuclei 
and n pairs can be written in the following form 



dr 



n+2 



n 

i=l 

n 

n 



d^q+d^q. 



fj^\^(27r)4 2 'J2si27rr 



(29) 



with q = p[j_. 

Cross section of n pairs production has the form 



dan 



1 \Mn\' 
8s n\ 



n+2- 



(30) 
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Statistical factor l/n\ is included to take into account the identity of pairs. 
Using the expressions (^)- (|29|) we get 

§!^ = P„(p), P„(p) = (^-^^(^))" e-^^.(p), n>2 (31) 

with 



Mp) = \B{p,r)\'ds^d'vd^r. (32) 

It can be easily recognized that for A{p) from (^) 

A(p) = yli(p), (33) 

thus confirming the Poisson character of probability distribution in impact- 
parameter representation. 

Let us mention that the effect of eikonal photons as well as the imaginary 
part of the amplitude corresponding to LBL blocks do not modify the total 
cross section as well as differential cross section integrated over phase volume 
of final nuclei. Really, integrating the square of the amplitude (^Bj) over the 
phase volume one immediately obtains 

3iq(Pi-P2)ei{[^(pi)-V'(P2)]-i[v(Pi)-¥'(P2)]}j(^p^)j(p2) — ^^^^ = 

TT TT TT 

/^(P)^- (34) 
7r 



Nevertheless, the exclusive cross section is sensitive to both these factors. 

Expression (|3^ can be simplified if one neglects the dependence of Comp- 
ton tensor ( [T^ ) on external photons virtualities 

B^p{k,r)e''e^ ^ 5„^(0, r)e"e'^ (35) 

and use the well known relation 

oo 

/EK(0,OeVf^*^3^. ,3a) 

Am? 

As a result the expression (^) can be cast in the form 

A{p) = Y^2i^'Z,Z,Mp), (37) 
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|ki| |ki — r| |k2| |k2 — r| n n n 

|k| |k — xl |k'| |k' + xl TT TT TT 

e^xPin2 ^ ^ 1^ (in(p^) + 0(1)) , 

V X / ^ P 

where we introduce the cut-off parameter < m as a result of the fast 
decreasing of matrix element of pair production by two photons. For the 
case of heavy leptons production (/i or r) the upper hmit must be replaced 
by quantity Q which can be associated with maximal momentum transferred 
to nucleus without its disintegration. For the case pm 3> 1 one has 

e;« (^2 y y Vl 

A{p) ^ (In(pm) + 0(1)) (38) 

which is in the agreement with [Q. 

Our formula (|3l| ) can be applied for the case n = (the probability of 
elastic nuclei scattering). The case n = 1 needs a bit accurate consideration. 
The expression for cxi can be cast in the following form 



a^ = aB + L j A{p) (e'^^^'') - l) d^p. 



(39) 



First term ob corresponds to the leading order of the Racah formula (see 
for instance [Q) for the cross section in Born approximation. Inferring (PUj), 
one has to take into account the longitudinal components of momenta of 
exchanged photons which create the pair. The second term takes into account 
the unitarity corrections to the total cross section. 



Conclusions 

In conclusions let us enumerate very shortly the main results of the present 
work. In our paper we obtained the general form for the amplitude of n 
lepton pairs production, accounting for the mutual Coulomb interaction of 
relativistic ions. We get the probability of vacuum-vacuum transition in the 
closed analytic form. We confirmed that the probability of n pair production 
can be approximated by Poisson distribution. Although this result has been 
claimed before in a number of publications, the mutual interaction of ions 
has not been taken into account in these works. At last we have shown that 
the Coulomb interaction between the colliding nuclei, although important in 
the exclusive kinematics, does not affect the integral yield of lepton pairs. 
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